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Abstract 

Depending on various assumptions on the energy scale of infla- 
tion and assuming a primordial power spectrum of a step-like 
structure, we explore new possibilities for Primordial Black Holes 
(PBH) and Planck relics to contribute substantially to Cold Dark 
Matter in the Universe. A recently proposed possibility to pro- 
duce Planck relics in four-dimensional string gravity is consid- 
ered in this framework. Possible experimental detection of PBHs 
through gravitational waves is also explored. We stress that infla- 
tion with a low energy scale, and also possibly when Planck relics 
are produced, leads unavoidably to relics originating from PBHs 
that are not effectively classical during their formation, rendering 
the usual formalism inadequate for them. 
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1 Introduction 



The formation of Primordial Black Holes in the early stages of the universe 
is a generic feature and it is therefore interesting to study its cosmological 
consequences pQ . Whatever the formation model, the Primordial Black Hole 
(PBH) spectrum must be in agreement with two types of constraints. The 
first one is associated with evaporation: the density must be low enough 
so that physical effects due to the Hawking radiation do not contradict any 
observed phenomena. They are based on the entropy per baryon, the nn pro- 
duction at nucleosynthesis, the deuterium destruction, the Helium-4 spalla- 
tion j2] [3] and, finally, on the nowadays observed gamma-ray [3] and antipro- 
ton |5] spectra. Those constraints apply for initial PBH mass between 10 9 g 
and 10 15 g, as these are the initial masses that can influence the above obser- 
vations through their (Hawking) evaporation. Using the quantity (3, which 
gives the probability that a region has the required density contrast to form a 
PBH at the horizon crossing time corresponding to the considered scale, the 
cosmic-ray constraints are the more stringent ones, leading approximately to 
P(Mpbh = 5 x 10 14 g) < 10~ 26 The second type of constraints is associ- 
ated with the normalization of the spectrum on cosmological scales probed 
by Cosmic Microwave Background (CMB) data. Whatever the considered 
power spectrum to form PBHs, it must generate a correct density contrast on 
COBE scales, i.e. large angular scales corresponding to the present Hubble 
radius scale. 

A pure scale- invariant Harrison-Zel'dovich power spectrum from the CMB 
scales up to very small scales would lead to a negligible amount of PBHs jB] , 
The only way to produce PBHs gnificant dark matter candidate is 

to increase the power on small scales without contradicting the observational 
data. A first attempt in this direction would be to allow for a tilt: P(k) oc k n 
with n > 1. Even without considering possible inconsistencies with cosmic- 
ray data, the required value, around n ~ 1.3 |2J, seems extremely disfavoured 
by the analysis of the most recent CMB experiments: between n « 0.91±0.06 
jH] (WMAP measurements: CMB + running spectral index) and n « 1.04 ± 
0.12 jU] (Archeops measurements: CMB + H Q ). A natural alternative is to 
boost power on small scales by means of a bump in the fluctuations power 
spectrum as suggested, for example, in [7j, [TU]. We follow here this idea and 
show that a (very) wide new parameter space can be opened for dark matter 
when the energy scale of inflation is low enough. Furthermore, in contrast 
with previous models, in such a scenario a simple step-like structure in the 
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spectrum is enough to generate a large quantity of PBHs without being in 
conflict with observations. Indeed, as a sufficient low energy scale allows 
the model to evade the 7-ray background constraints, no bump is required 
in the mass variance. Some observational probes through the emission of 
gravitational waves by coalescing PBHs are suggested to test this hypothesis. 
The paper deals also with new models for Planck relics formation, based on a 
four- dimensional effective action in the framework of string gravity, to show 
that even for a high-energy scale inflation, PBH-induced dark matter is a 
viable candidate. We stress that in some of these scenarios the production 
of (non-evaporated) PBHs from quantum fluctuations which are not highly 
squeezed and therefore not effectively classical [TT] . is unavoidable and cannot 
be handled with the usual formalism ^2]- This will be the case when the 
energy scale of inflation is low enough so that the Hubble mass at the end of 
inflation, Mn, e , is larger than M* « 5x 10 14 g, the mass of PBHs ending their 
evaporation at the present time, and also in high-scale inflation for Planck 
relics whose initial mass at formation is close to M# e <C M*. 

2 Inflation with a low energy scale 
2.1 A pure step 

An important consequence of low scale inflation is the decrease of the re- 
heating temperature. Though in practice the reheating scale can be much 
lower than the energy scale at the end of inflation, it will be enough for 
our purposes to make the simplifying assumption that the reheating is in- 
stantaneous. A low reheating temperature is required in order to avoid the 
possible overproduction of gravitinos 13 : Trh < 10 8 GeV. This value was 
even decreased to 4 x 10 6 GeV in some works based on Lithium abundance 
|14j . This makes the horizon size at the end of inflation very large with an 
associated Hubble mass Mu > 10 16 g. This point is extremely important for 
PBH dark matter as it allows to avoid the main problem explained in e.g. 
|15j . namely the gamma-ray constraint which comes from the contribution 
to the 7-ray background of evaporated PBHs. Then only the gravitational 
constraint, namely constraints on the present abundance of PBHs, would 
apply for PBH masses Mpbh greater than M^ e , the Hubble mass at the end 
of inflation. 

Previous work on the subject argued that one way to produce a significant 
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amount of dark matter in the form of PBHs is to increase the mass variance 
in a well localized region so as to remain in agreement with the gamma-ray 
constraint. For example, using the Broken Scale Invariance (BSI) Starobin- 
sky spectrum [T^], it was shown that the oscillation in the power spectrum 
due to the jump in the derivative of the inflaton potential should produce a 
bump in the mass variance ^3] ■ This slight increase in variance o 2 H can boost 
the PBH formation probability (3 by more than ten orders of magnitude. The 
resulting bump in the probability (3 to form PBHs can yield flpBH,o — 0.3 
for 5 x 10 15 g < Mpbh < 10 21 g with values of p « 8 x 10~ 4 [15J where p 2 is 
the ratio of the power on large scales with respect to that on small scales. 

Clearly, if the horizon mass Mji, e & t the end of inflation is larger than 
M* ~ 5 x 10 14 g, the initial mass of a PBH whose lifetime if equal to the age 
of the Universe, the 7-ray and antiproton constraints as well as all the other 
constraints on smaller masses associated with evaporation are automatically 
evaded without any requirement about the shape of the fluctuation spectrum. 
PBHs with masses above M* are nearly insensitive to the Hawking emission 
as the temperature T = he 3 / {^GkMpBu) becomes smaller than the rest 
mass of any known massive field. An extremely wide mass range without 
constraint (except, to some extent, for microlensing upper limits) is therefore 
opened. 

The relative PBH abundance today is given by [7j 

ttpBH. (MpBH) « 1-3 x 10 17 (3(M PBH ) , (1) 

\MpBH J 

for h ~ 0.7, and the subscript stands for the present-day value. The 
quantity (3 is defined as 



'max S 



P(M H ) = ~ 1 = 1 / e 2 W d5 sa J^ k) e , (2) 



where is the horizon crossing time for the considered mode, S is the density 
contrast, M H is the Hubble mass at % and cr 2 H {t k ) = u 2 (R) \t k - For an accurate 
calculation, the mass variance & 2 (R) =< ($rf)% > is computed with a top- 
hat window function Wth and R = ^—\t k , using the expression [T%] 

fee 

°H{tk) = i / k ^ F{kx) W 2 H (c s x) W 2 H {x) dx . (3) 

0-L 71 " Jo 



3 



In (J3J), k e corresponds to the Hubble crossing scale at the end of inflation, 
c 2 s = |, the quantity F{k) is defined through the equality k 3 P(k,tk) = 

(|) 4 F(A;) for tk < t eq where P(k,t) is the power spectrum of the primor- 
dial fluctuations [H],[IH]- However, in order to give a conservative estimate 
of the increase in power on small scales, we can assume following [7j, that a 
(possibly smoothed-out) jump occurs around some characteristic mass M s , 
Mn, e < M s -C MH{t eq ), in the mass variance spectrum, viz. 

a c H OBE = p a H {t k ) . (4) 

Normalizing the mass variance with the CMB large angular scales measure- 
ments, the increase in power and therefore the amplitude p of the step can 
be estimated as a function of the horizon mass at the end of inflation M# e 
as follows 



V 



^COBE 
a H 



5, 



mm \ 



LW 



1.7 x 10 34 



2nQ 



PBHfi L lvl H,e 



1015 g 



(5) 



where LW stands for the Lambert- W function (with LW(xe x ) = x). With 
QpBH,o ~ ^m,o ~ 0.3, the numerical estimates are : p ~ 6.5 x 10~ 4 for 
M Hj£ = 10 15 g, p ss 5.5 x 10~ 4 for M H)E = 10 25 g, p ss 4.1 x 10~ 4 for 
M H e = 10 35 g. In these estimates we have taken 5 min = 0.7. In principle, the 
reheating temperature can be as low as the MeV scale (the nucleosynthesis 
temperature), leading to huge horizon masses around 10 38 g. This can be 
considered as the upper limit for the low-mass cutoff of PBH spectra. It is 
interesting that this corresponds to the highest viable PBH masses if CDM 
is made of PBHs [T^. It opens a very wide parameter space (M# )e ,M s ) for 
PBH dark matter. Furthermore, if the reheating temperature is smaller than 
1 GeV (M HjB > 10 16 g), PBHs could be one of the viable CDM candidates 
left, as super symmetric dark matter cannot contribute substantially to dark 
matter 



2.2 Gravitational waves 

Probing PBHs as a dark matter candidate is experimentally is very difficult. 
As long as their masses are greater than 10 15 g, black holes do not radiate and 
become really black. A decisive way to detect them, and to observationally 
confirm or exclude this model, could be to look for gravitational waves from 
coalescing PBHs. The maximum distance R ma x between the Earth and the 
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binary system compatible with the sensitivity of a given detector for a fixed 
PBH mass M PBH is given by [2$] : 
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20 Mpc J 



where hsBmin is the minimum characteristic amplitude of a wave lying within 
the detectot sensitivity and v is the considered frequency. To estimate the 
number n(M PBH , R max ) of PBHs inside this sphere, we used an isothermal 
profile inside the Milky- Way halo (for R < 150 kpc) which is given by : 



R c + R & 
' R 2 C + R 2 Q - 2rR Q cosij + r 2 



where p ~ 5 x 10~ 25 gcm~ 3 is the local halo mass density, Rc ~ 3 kpc is the 
core radius, R® ~ 8 kpc is the galactocentric distance, r is the distance of 
the binary system to the Earth, ip is the angle between the considered point 
and the galactic center seen from the Earth, and finally p(r, iji) is the halo 
mass density at coordinates r, ip. This leads to: 

„ v irp Q R 2 c + Rl [ Rma \ ( R 2 c + R 2 Q + 2rR Q + r 2 ) , 
n{M PBH , R max ) = — ^ C D ° / In <^ ^ ° ~ ^ rrfr. 



Mpbh R G Jo \R 2 c + R%-2rR Q + r 2 

(8) 

For ^> 150 kpc, an average dark matter distribution with p m 0.3p c is 
assumed. Finally the coalescence rate / within this volume is computed 
under the natural assumption that the distribution function of the PBHs 
comoving separation is uniform and can be straitfowardly obtained from |3*U] : 

f „ 3 ( m pbh \ * n(M PBH} R max ) 
V M J t 

where to is the age of the Universe. Gathering all those formulae and perform- 
ing numerical estimates shows that if PBHs have masses above 2 x 10~ 5 M Q , 
they should generate more than one "event" per year in the VIRGO detec- 
tor. If the LISA frequencies and sensitivity are considered, the minimal mass 
decreases down to 1O _11 M . The interesting mass range probed covers then 
nearly fifteen orders of magnitude, though it also overlaps with microlensing 
data which exclude a significant contribution between 2 x 1O~ 7 M and 1M 
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2.3 "Quantum" relics 

A generic feature of low scale inflationary models is that all PBHs that are 
produced will survive and not evaporate. In particular, those PBHs that 
would form right after inflation correspond to fluctuations that were not (long 
enough) outside the Hubble radius and are therefore not highly squeezed. 
Indeed, for the first scales which re-enter the Hubble radius, producing in 
these models PBHs with Mpbh > Mn,e, the fluctuations are not highly 
squeezed. This is equivalent to saying that the decaying mode is still present 
and actually of the same order as the growing mode, and cannot be neglected. 
Hence such low scale inflation models lead to the possible production of 
PBHs by inflationary fluctuations that cannot be considered as stochastic 
classical fluctuations, so the PBHs are not evaporated by today, in contrast 
to high scale inflation where Mn, e *C M* and for which PBHs with mass 
Mpbh ~ Mn t e have evaporated already during the (very) primordial stage 
of our Universe. Even if there is no significant increase in power on small 
scales, PBHs will be produced that cannot be described as classical objects 
and that will survive till the present time. We might therefore call them 
quantum relics. The intriguing point is not their abundance, which should 
be low, but rather the very nature of these surviving objects. 

We want to illustrate the production of "quantum relics" with a concrete 
high energy physics inspired low scale inflationary model. This production 
can take place, for example, in the phenemenological model considered in the 
first part of this Section. 

It would be interesting if the scale of inflation corresponds to the super- 
symmetry breaking scale or even the electroweak scale. Initial conditions 
(through thermal effects) could set the inflaton field <ft close to the origin 
where some symmetry is unbroken, as in "new inflation". Inflation then 
takes place at small field values. At low temperature the inflaton starts 
rolling away from the origin = (the effective mass term becomes nega- 
tive), spontaneously breaking the underlying symmetry. The following quite 
general inflationary potential can be considered in the context of supergravity 
inflation |2T] : 



where A is the near-constant vacuum energy driving inflation, M p the Planck 
mass, b a "bare" mass term, c a logarithmically varying mass term brought 




(10) 
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by radiative corrections, and k, p and q determine the end of inflation by 
inclusion of higher-order terms. It has been shown that a value for A from 
as low as 1 GeV up to ~ 10 11 GeV can be obtained for reasonable choices of 
p and q, while still generating an acceptable spectrum of perturbations with 
a spectral index n s < 0.95 close to 1 and a sufficient number of e-folds. So 
the Hubble mass M^ e at the end of inflation will satisfy M^ e > M*. We 
can quantify the degree of classicality of the formed PBHs with the ratio D 
of the growing to the decaying mode for the scales under consideration [T2*] . 
A ratio D ^> 1 for a given scale corresponds to effective classicality of the 
fluctuations on this scale. We estimate it for adiabatic fluctuations in these 
models and we find towards the end of inflation, following 12 j, the expression 



where M 4 stands for the energy density at the time when the PBH is formed. 
It is related to Mh through 



The expression (jllj) assumes that the inflationary energy density scale ~ A 4 , 
and it applies to PBHs formed shortly after the inflationary phase. As can be 
seen from pijl . only those PBHs produced immediately after the inflationary 
phase have D ~ 1 and are therefore quantum objects. Though PBHs are here 
clearly irrelevant as CDM candidates because of their negligible abundance, 
the appearance of these "quantum relics" might be one of the few ways in 
which the quantum nature of the inflationary fluctuations can be exhibited. 

3 Inflation with a high energy scale 
3.1 Two distinct inflationary stages 

Another way to evade the "small scales" problems for PBH formation still 
in the framework of high scale inflation is through the existence of a second 
inflationary stage at much lower energies. So, a first stage of inflation solves 
all the problems usually solved by inflation, generates the cosmological per- 
turbations observed, and produces also PBHs including in the "dangerous" 




(11) 




(12) 
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mass interval around M*. However, due to the second stage of inflation, a 
significant PBH abundance produced during the first inflation is allowed. We 
give now in full generality the salient features of such a scenario. 

Let us assume that the second inflation starts when the Hubble mass 
equals Mha, at a much lower energy than the first inflation scale. For sim- 
plicity we can assume the Hubble mass is constant during the second inflation. 
On a large range of scales, fluctuations produced during the first inflation will 
re-enter the Hubble radius, thereby possibly producing PBHs. Part of those 
scales which re-entered the Hubble radius between the two inflationary stages 
will be expelled again outside the Hubble radius during the second inflation. 
Let us consider the scale k H = aiHi that corresponds to the Hubble radius 
at the beginning of the second inflation. It will eventually, at the time tj, H , 
reenter the Hubble radius when the Hubble mass is given by Mn{tk H )- It is 
easy to derive the following relation between MH{tk H ) and Mha, 

M H (t kH ) 2N 



M, 



(13) 



HA 



where N is the number of e-folds during the second inflation. Due to a 
much lower energy scale, the amplitude of the produced fluctuations is quite 
negligible and will not lead to a significant PBH production. Only the fluc- 
tuations of the first, high scale, inflation will. Therefore there will be a gap 
in the mass range Mha < Mpbh < Mnitku). In addition, the density of all 
objects created before the second inflation will be reduced by an additional 
factor e~ 3N . Hence, the only significant abundance of PBHs corresponds to 
the range 

M PBH > e 2N M H:i . (14) 

Clearly it is possible to have PBHs as CDM in this range and still evade the 
small scales constraint coming from the evaporated PBHs. 

An interesting low scale inflationary model of that kind is thermal infla- 
tion [22] , triggered by a scalar field termed flaton, which can appear in SUSY 
theories. The consequences of thermal inflation on PBH abundance were con- 
sidered in |231 • By definition, a flaton has a large vacuum expectation value 
M > 10 3 GeV, while having a mass of order the electroweak breaking scale 
m ~ 10 2 — 10 3 GeV. This leads to an almost flat potential for the flaton field 
/: V ~ Vo — m 2 \f\ 2 with Vo — m 2 M 2 . During thermal inflation the flaton 
field is held at the origin by finite temperature effects and the potential is 
dominated by the false vacuum energy Vq. Thermal inflation starts at the 



8 



temperature Tj ~ \JmM when the thermal energy density falls below Vo and 
ends at T e ~ m when the flaton can escape the false vacuum. The number 
of e-folds N is then immediately given as N — | In® and the density of 
all PBHs produced before thermal inflation is suppressed accordingly by a 
factor (§)§. If we take M ~ 10 11 GeV, m ~ 10 3 GeV, thermal inflation 
starts at Tj ~ 10 7 GeV which corresponds to a Hubble mass Mh,% ~ 10 18 g, 
and ends at T e ~ 10 3 GeV. The number of e-folds is iV « 9 and the density 
of PBHs with mass M < M B ,i is suppressed by a factor ~ 1CT 12 . Clearly it is 
then possible to have a significant amount of PBHs and still evade the small 
scales constraints. Note however that with the numbers given above, PBHs 
as CDM can only exist in the range Mpbh ^ 10 26 g, starting near the edge of 
the range probed by the EROS data which constrain galactic dark matter in 
the range 2 x i(T 7 M < Mpbh < 1M |2H- In particular, this could apply to 
the model with a step considered in the previous chapter, followed by a stage 
of thermal inflation. However, a step in the primordial spectrum produced 
during the first period of inflation at the characteristic scale M s < 2 x 1O~ 7 M 
leaves possibly only a tiny mass interval (M H (tk H ) < M PBH <2x 1CT 7 Mq) 
where PBHs are not constrained by observations. There will be essentially 
no "quantum relics" in this scenario as practically no PBHs are produced 
during the second (low-scale) inflation. 



3.2 Planck relics 

If we simply ignore moduli and gravitinos that appear in supersymmetric 
theories and their possible overproduction if the reheating temperature is 
too high, as none of them as ever been detected, still another interesting 
possibility could be the production of Planck relics. Indeed, the unavoidable 
upper limit is imposed by gravitational waves and fixes the smallest possible 
horizon mass after inflation around 1 g. If the horizon mass is in this range, 
a natural way to produce dark matter is through PBH relics. The idea was 
first mentioned in |23|. Nevertheless, two critical ingredients were missing 
at that time : the normalization of the primordial spectrum to COBE data 
and a realistic (or, at least, possible) model to stop Hawking evaporation in 
the Planck era. The latter point received a new light in the framework of 
string gravity. The following four- dimensional effective action with second 
order curvature corrections can be found: 




-R + 2d^<f) + \e- 2 *S GB + 



(15) 
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where A is the string coupling constant, R is the Ricci scalar, is the dilatonic 
field and Sqb — RijkiR 1 ^ 1 ~ ^RijR 1 -' + R 2 is the curvature invariant - the 
Gauss-Bonnet term. This generalisation of the Einstein Lagrangian leads to 
the very important result that there is a minimal relic mass M re \ for black 
holes produced in this model • Solving the equations at the first order of 
perturbation gives a minimal radius: 



rT f = y/X V 4 ^(^°), (16) 

where <fih{<fioo) is the dilatonic value at the horizon. The crucial point is that 
this result remains true when higher order corrections or time perturbations 
are taken into account. The resulting value should be around 2M P . It is 
even increased to 10M p if moduli fields are considered, making the conclu- 
sion very robust and conservative. The subsequent decrease of the Hawking 
evaporation leads to an asymptotically stable state |2Zj, giving a quantitative 
argument in favour of the existence of Planck relics. For formation masses 
above 10 9 g, important constraints are associated with Helium and Deu- 
terium destruction. Actually, as we will see below, the relevant upper limit 
for the initial PBH mass is ~ 10 5 g. 

Once again, we consider a mass variance spectrum with a characteristic 
scale M s . For small initial masses Mpbh ~ 1 g, the Planck relics relative 
density Q re i can be written as, cf.(l), 



M rel /10 i& g 



n^o « 1-3 x 10^(3(M PBH )-^- ( ) (17) 

MpBH \MpBH J 

3 

/ 1 f)15 „ \ 2 

« 2.83 7 x 10- 3 /3(ilW) — , (18) 

\MpBH J 

where we have taken M re \ = 7M P , while Mpbh refers to the initial PBH 
mass. This leads to the following value for the step amplitude: 



V 



.COBE 

u 



® H 



A 



LW 



8.0 x 10" 6 



M 



rel 



1015 g 



(19) 



If Planck relics are to explain fl mi o ~ 0.3, p varies from 7.1 x 10~ 4 to 5.5 x 10~ 4 
for initial PBH masses between 1 g and 10 5 g with M re i = M p and between 
7.3 x 10 -4 and 5.6 x 10 -4 with M rd = 10M„. As for non-evaporating PBHs, 
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p has a very slight dependence on Mjj, e as (3 is extremely sensitive to the 
mass variance in this range. 

Equation (JT8J) is very accurate for small masses Mpbh - 1 g and its 
validity extends up to M PBH ~ M' = 8 75 x 10 5 g, with the corresponding 
range 1CT 21 7 _1 < (3{Mpbh) % 10 -12 7 _1 . We have checked this numerically 
using the simplifying assumption that all the evaporation produces either a 
relativistic, or else a non-relativistic, component. For Mpbh up to M', the 
energy density due to evaporation is still much smaller than the preexisting 
radiation background, at the time when the relics have formed. However, it 
should be stressed again that when Mpbh ~ Mn, e , the quantity (3 entering 
(fTHj) loses its meaning as a probability, not to mention the fact that the 
asymptotic domination of the growing mode is not achieved in this regime. 
Beyond M', PBHs dominate the energy density before their evaporation is 
completed. This gives rise to a different expression for fl re j which is nearly 
independent of (3 (.3, 28J). Interestingly, in this mass range it is possible 
to obtain PBHs as CDM in significant amounts only for masses Mpbh — 
M'. Thus, in the context of CDM, only equation (jTHj) is relevant and PBHs 
cannot contribute significantly to CDM if their mass lies in the range M' < 
Mpbh < 10 9 g. Surprisingly, in such a scenario too, the notion of quantum 
relics resurfaces for those Planck relics originating from initial PBH masses 
Mpbh — Mn, e ~ 1 g. Because of their supposed large abundance, Planck 
relics with initial mass Mpbh ~ Mn, e ~ 1 g would probably be ruled out 
but we conjecture that a bump producing PBHs around 10 5 g could yield a 
viable CDM candidate. 



4 Conclusion 

In this work we have presented results of two different kinds. The first 
kind refers to the possibility of having PBHs as a viable CDM candidate. 
We have shown that PBHs can be produced in significant amounts when 
the inflationary energy scale is low enough without being in conflict with 
constraints coming from the evaporation of PBHs with masses M < M* 

5 x 10 14 g, simply because in these models such PBHs do not form. On the 
other hand, if the inflationary scale is high, PBHs can still remain viable CDM 
candidates in theories where Planck relics exist, and we have considered this 
possibility too. Still another possibility is to have two inflationary stages, 
a high energy scale inflation followed by another inflation on a much lower 
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energy scale. The possibility to detect PBHs through gravitational waves 
was also considered. 

The second kind refers to the unavoidable production in low energy scale 
inflation of PBHs that form from genuinely quantum-mechanical fluctuations 
and that are not evaporated by today. We have called these objects "quantum 
relics". As we have pointed out, such objects can also form in the context of 
Planck relics. Clearly, the very existence of such objects poses a number of 
fascinating open problems both with respect to their formation and to their 
subsequent evolution and possible decoherence (see e.g. [HI]). 
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